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Abstract 

We analyze the theoretical uncertainties in Br(B — > X s 7) due to the choice of the high 
energy matching scale nw = C(M W ) and the scale fi t at which the running top quark mass is 
defined: m t (fi t ). To this end we have repeated the calculation of the initial conditions confirming 
the final results of Adel and Yao and Greub and Hurth and generalizing them to include the 
dependences on \i t and fiw with \i t ^ fi w . In the leading order the fj,yy and fi t uncertainties 
in Br(B — ► X^) turn out to be ±13% and ±3% respectively. We show analytically how these 
uncertainties are reduced after including next-to-leading QCD corrections. They amount to 
±1.1% and ±0.4% respectively. Reanalyzing the uncertainties due to the scale = 0{m b ) we 
find that after the inclusion of NLO effects they amount to ±4.3% which is a factor 2/3 smaller 
than claimed in the literature. Including the uncertainties due to input parameters as well as the 
non-perturbative l/m| and 1/m^ corrections we find Br{B^X s r )) = (3.60 ± 0.33) x 10 -4 where 
the error is dominated by uncertainties in the input parameters. This should be compared with 
(3.28 ±0.33) x 10~ 4 found by Chetyrkin et al. where the error is shared evenly between the scale 
and parametric uncertainties. 
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1. The inclusive B — > X s 7 decay has been subject of considerable experimental and theoretical 
interest during the last ten years. Experimentally its branching ratio is found by the CLEO 
collaboration to be pj 

Br(B -> X sl ) = (2.32 ± 0.57 ± 0.35) x 10" 4 , (1) 

and a very recent preliminary result from the ALEPH collaboration reads 

Br(B -> X s7 ) = (3.38 ± 0.74 ± 0.85) x 10" 4 . (2) 

In ([!]) and (§) the first error is statistical and the second is systematic. On the other hand the 
complete NLO analysis gives 

Br(B^X sl ) = (3.28 ± 0.22 (scale) ± 0.25 (par)) x 10~ 4 = (3.28 ± 0.33) x 10~ 4 . (3) 

where the first error results from the scale uncertainty (see below) and the second error from the 
uncertainties in the input parameters. A similar result has been obtained in 0j. 
The NLO analyses presented in || £| reduced by a factor of 3-4 the /^-uncertainties || |J present 
in the leading order, where //& = 0(m&) is the scale at which the relevant decay matrix element 
is evaluated. This reduction of the /!& uncertainty is very welcome because in the forthcoming 
years much more precise measurements of Br{B^-X s r y) are expected from the upgraded CLEO 
detector, as well as from the B-factories at SLAC and KEK. This is also the reason why continuing 
efforts are being made to estimate non-perturbative corrections to the B—^-Xsj decay with higher 
precision [8], || |10[ as well. It appears that these latter corrections amount only to a few percent 



and constitute a rather small theoretical uncertainty. 

2. In this letter we have repeated the numerical analysis of || to find that the remaining scale 
uncertainties are by roughly a factor 1.5-2.0 smaller than quoted by these authors and in [|J. 
This includes also two additional theoretical uncertainties which have not been addressed in 
the literature. They are related to the choice of the high energy matching scale fi\y = C(M W ) 
and the scale \xt = 0(m t ) at which the running top quark mass is defined: m t (/i 4 ). These two 
scales enter the analysis of B—^Xsj in the process of calculating the initial conditions for the 
renormalization group running of the Wilson coefficients Cj and C% of the operators 

= ^m b s a a^(l + ^)b a F^ Q 8 = ^m b s a a^(l + l5 )T^b p G%. (4) 

Here e and g s denote the electromagnetic and strong coupling constants respectively. These initial 
conditions have been calculated at NLO in [|TTJ and have been recently confirmed in These 



NLO corrections are necessary to remove the renormalization scheme dependence present in the 
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renormalization group evolution from fi w = C(M W ) down to fif, = 0(m h ). From our point of 
view the additional reason for performing these rather tedious calculations is the reduction of the 
uncertainties related to the choices of and \i t . These uncertainties have not been discussed 
in [0, 



To this end we have repeated the calculation of the initial condition for the by far dominant 



Wilson coefficient C-j confirming the final result in flTij , |12|| and generalizing it to include the 
dependences on \i t and nw with \i t ^ fiw In 101 anc ^ 13 ^ w = A 4 * = Vwt have been used. The 
technical details of our calculation which differs in certain aspects from the previous ones will be 



presented elsewhere |L3] . Here we discuss first the issue of the nw and /i t uncertainties and their 
reduction after the inclusion of NLO corrections. Subsequently we discuss the //& uncertainties 
and we present our estimate of BriB — > X s ^) in the Standard Model. 
3. In the leading logarithmic approximation one has 

Cf^l 2 , (5) 
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m b,pole 



is the phase space factor in Br(B — > X c ev e ) and a = e 2 /47r. 

The effective renormalization scheme independent coefficient Cy (^t&) introduced in || is given 
by 

C^M = V*C?\» W ) + ? (rjM - V M) C f\» w ) + £ , (7) 
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with 



(o), 3x^ - 2x 2 -8x 3 - 5x 2 + 7x t , * 

° 7 {fIw) ~ A(x t - IY ln Xt + 24(x t - 1)3 ' (9) 

^(o), x -3x t 2 -x 3 t +5x 2 + 2x t . . 

C * M = M^W lUXt+ 8(3,-1)3 (10) 

Xt = (11) 
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The numbers a, and hi are given in table [I]. 
There are three scale uncertainties present in fl5|): 



The low energy scale /if, = 0(m b ) at which the Wilson Coefficient (fJ-b) is evaluated. 
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Table 1: Magic Numbers. 

• The high energy scale fiw — 0{M^) at which the full theory is matched with the effective 
five-quark theory. In LO this scale enters only r/. C^il^w) an d Cf\nw)i usually denoted 
by cf^Mw) and Cf } (M W ), serve in LO as initial conditions to the renormalization group 
evolution from fiw down to fib- As seen explicitly in (H) and ([U]) they do not depend on 

• The scale = 0(m t ) at which the running top quark mass is defined. In LO it enters only 
x t in ©. 

It should be stressed that fiw and /i f do not have to be equal. Initially when the top quark and 
the W-boson are integrated out, it is convenient in the process of matching to keep fi t = Hw- 
Yet one has always the freedom to redefine the top quark mass and to work with m t (fi t ) where 
7^ fJ'W- It is evident from the formulae above that in LO the variations of fib, H>w an d 
fit remain uncompensated which results in potential theoretical uncertainties in the predicted 
branching ratio. 

In the context of phenomenological analyses of B — >• Xg'-f, only the uncertainty due to fib has 
been discussed ]5], || ^ f| . It is the purpose of this letter to analyze the uncertainties due to fiw 
and fi t and to reanalyze the /ib-uncertainty. 

It is customary to estimate the uncertainties due to fib by varying it in the range mb/2 < fib < 
2mb. Similarly one can vary fi\y and fit in the ranges Mw/2 < fiw < 2M\y and m t /2 < fi t < 2mt 
respectively. Specifically in our numerical analysis of fiw and fi t uncertainties we will consider 
the ranges 

40 GeV < fi w < 160 GeV 80 GeV < fi t < 320 GeV (12) 
setting fi b = nib = mb, P oie = 4.8 GeV. 
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In the LO analysis we use 



with a s (M z ) = 0.118 and 



v(n) 



a s {M z ) (M z 
v(ii) = 1 - (Jo — In 



2tt 



m t (fi t ) = m t (m t ) 



(13) 



(14) 



a s (m t ) 

We are using the parameters and rh^ defined in the effective theory with five flavours 
throughout in this work, hence j3 = 23/3. We set m t (m t ) = 168 GeV and m t = m tyP oie = 
176 GeV. 

Varying and fi t in the ranges ([TJ) we find the following uncertainties in the branching ratio: 



ABr(B -> V s7 ) 



(15) 



±13% (n w ) 
±3% (//*) 

to be compared with the ±22% uncertainty due to the variation of the scale //& 0, 0]. The fact 
that the /ziy-uncertainty is smaller than the /i& uncertainty is entirely due to a s (/i^) < ot s (ii b ). 
Still this uncertainty is rather disturbing as it introduces an error of approximately ±0.40 • 10~ 4 
in the branching ratio. The smallness of the /^-uncertainty is related to the weak x t dependence 
of Cf\jj,w) and Cf\jj,w) which in the range of interest can be well approximated by 



-0.122 x 



0.30 



-0.072 x 



0.19 



(16) 



Thus even if 161 GeV < m t (/z t ) — 178 GeV for fx t in (|i2|) , the \x t uncertainty in Br{B — > X s ~f) 
is small. This should be contrasted with B s — > fip,, Kl is^vv and B° — B° mixing, where \it 



uncertainties in LO have been found [FJ, |Tj|] to be ±13%, ±10% and ±9% respectively. 

4. We will next investigate how much the uncertainties in ( 115] ) are reduced after including NLO 

corrections. 

The formula @ modifies after the inclusion of NLO corrections as follows ||: 
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k(z) \ m feiPole J k(z) 
with k(z) being the QCD correction to the semileptonic decay |T6[ and given to a good approxi 
mation by fli"7 



k(z) 
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(19) 



5 



An exact analytic formula for k(z) can be found in [18]. Here Jib = 0{m^) is a scale in the 
calculation of QCD corrections to the semi-leptonic rate which is generally different from the one 
used in the b — > 57 transition. In this respect we differ from Greub et al. M who set /2& = fib- In 
[Q the choice fib = wfy has been made. We will return to this point below. 
Next 



D = Cf^^b) + 



47T 



i=l 



(O)eff 2 
, Ti7 1 m b 



(20) 



where C^ eS (^b) is the NLO correction to the effective Wilson coefficient of Qj\ 
Generalizing the formula (21) of ][3| to include ft, t and [l w dependences we find 1 
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where in the MS scheme 



Ci lW ^ w ) = C^(M w )+8x t 
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1 We would like to thank the authors of [|T9| for pointing out the missing logarithmic term in the original 
version of this work. See also the discussion after equation 
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The formulae for Cj 1 ^^ (My/) given above and presented in || are obtained from the results 
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IT] , |12]1 by using the general formulae for the effective coefficient functions |J. The formula 



for Cj^ (My/) has been confirmed by us [13 



The numbers e^-l{ are given in Table [I]. We have confirmed these numbers as well as the numerical 
coefficients in ( p2"D using the anomalous dimension matrices in ||. Next the 7] in (|8|) should now 
be calculated using the NLO expression 



a s (n) 



ot s (M z 



1 - 



A ol s (M z ) \nvQi) 



(28) 



where v(fi) is given in fllTf) and /?i = 

The constants resulting from the calculations of NLO corrections to decay matrix elements 
HJ are collected in where also explicit formulae for C| with z = 1 — 6, 8 and the values 

of 7j7' leff can be found. It should be stressed that the basis of the operators with i = 1 — 6 used 
in [^J differs from the standard basis used in the literature |2~3]] . For the discussion below it will 



be useful to have [24 



7 <r 



(29) 



416 (0)eff = 70 
81 728 27 
which enter (|23|) and (^4|) respectively. 

Finally the term A in flTT| ) originates from the bremsstrahlung corrections and the necessary 
virtual corrections needed for the cancellation of the infrared divergences. These have been 
calculated in [|(], |2l| and are also considered in ||, f| in the context of the full analysis. The 
explicit formula for A, which we use in our numerical analysis, can be found in equation (32) of 

10- 

Setting \iy/ = n t = /j, wt , replacing j^ eS by its value in the NDR scheme 727 >NDR = 464/81 and 
adding all fii dependent terms in (|^) we recover the [iwt dependence of C^ NDR (fiwt) found in 
|i~2"f . Similarly replacing 72g )eff by 728 >NDR = 76/27 in (|4]) we recover the [iy/t dependence of 
Cg 1 ' )NDR ( / upi/ t ) given in []I2"|| . For /j, w = jj, t = M w the formulae above reduce to the ones given in 



2 In the replacement version of || several quantities entering the formula for A have been corrected. In this 
paper the updated values are used. We thank M. Neubert [Q, P. Gambino and M. Misiak for informing us about 
these modifications. Accordingly the numerical results of this work are changed slightly. 
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5. Before entering the numerical analysis let us demonstrate analytically that the fi t and fi\y 
dependences present in C^ eS (fib) are indeed cancelled at 0(a s ) by the explicit scale dependent 
terms in fl2"3]). The scale dependent terms in fl2"3| ) do not enter this cancellation at this order 
in a s in B —>■ X^. On the other hand they are responsible for the cancellation of the scale 
dependences in C^ eS (fib) relevant for the b — > s gluon transition. 
Expanding the three terms in ([!]) in a s and keeping the leading logarithms we find: 

rpCfW) = fl + ?T ln 4) CPW) (30) 
\ 4vr 3 fi l w J 

\ (,» - „i) CPM = In if-C<°V„) (31) 
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respectively. In (|32|) we have used Y,hi — 0. Inserting these expansions into (p0|), we observe that 
the /ivi/ dependences in (|3"0"D, (pi|) and (^) are precisely cancelled by the three explicit logarithms 
in ( p3|) involving respectively. Similarly one can convince oneself that the /i r dependence of 



C^ eS (/j,b) is cancelled at 0{a s ) by the In fif/M^ term in (|23| ) 



Interestingly the last logarithm in (p2|) does not contribute to any cancellation of the /ivy 
dependence at this order in a s due to the relation Yh=i (| e « + = which can be verified by 
using the Table [L| 

Clearly there remain small fi t and fi\y dependences in flip) which can only be reduced by going 
beyond the NLO approximation. They constitute the theoretical uncertainty which should be 
taken into account in estimating the error in the prediction for BriB — > X^). 
Using the well known two-loop generalization of ( |T4D and varying fiw and m in the ranges (|l2l ) we 
find that the respective uncertainties in the branching ratio after the inclusion of NLO corrections 
are negligible: 

ABr(B^X,, ) = { ±Ll% M (33) 
\ ±0.4% (,,,) 

6. We have next performed the NLO analysis of the /j& dependence. Varying fib in the range 
2.5 GeV < fi b < 10 GeV we find 

ABr{B -> X s7 ) = ±4.3% (/i 6 ) (34) 

This reduction of the /i^-uncertainty by roughly a factor of five relative to ±22% in LO is caused 
by the presence of the explicit logarithm In ml / fj% in (p0|) . We note that our result in ([34]) differs 
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from the //&- uncertainty of ±6.6% quoted in JJ. A discussion with the latter authors confirmed 
our result. 

Next we would like to comment on the uncertainty due to variation of fit in k(z) given in 
( |19D . In |Q the choice flf, = fib has been made. Yet in our opinion such a treatment is not 
really correct, since the scale fib in the semi-leptonic decay has nothing to do with the scale 
fib in the renormalization group evolution in the B — ► X s 7 decay. Varying Jib in the range 
2.5 GeV < fi b < 10 GeV we find 

ABr{B -> X s7 ) = ±1.7% (/2 6 ) (35) 

Since the //& and /2& uncertainties are uncorrelated we can add them in quadrature finding ±4.6% 
for the total scale uncertainty due to fib and fib- This is smaller by roughly 30% than the case 



in which fib = fib is used. The addition of the uncertainties in fi t and fiw in (p3D modifies this 
result slightly and the total scale uncertainty in BriB — > X s ^j) amounts then to 

ABr(B^X s j) = ±4.8% (scale) (36) 

which is roughly by a factor of 1.5 smaller than quoted in ||, [|. 

It should be stressed that this pure theoretical uncertainty related to the truncation of the 
perturbative series should be distinguished from parametric uncertainties related to a s , the quark 
masses etc. discussed below. 

In our numerical calculations we have included all corrections in the NLO approximation. 
To work consistently in this order, we have in particular expanded the various factors in ( |TTD in 
a s and discarded all NNLO terms of order a 2 s which resulted in the process of multiplication. 
This treatment is different from [^, M, where the a s corrections in ( |T8"D have not been expanded 
in the evaluation of (|17D and therefore some higher order corrections have been kept. Different 
scenarios of partly incorporating higher order corrections by expanding or not expanding various 
factors in ( |T7D affect the branching ratio by ABr(B — > Xg'y) rs ±3.0%. This number indicates 
that indeed the scale uncertainty in fl5E]) realistically estimates the magnitude of yet unknown 
higher order corrections. 

The remaining uncertainties are due to the values of the various input parameters. In order 
to obtain the final result for the branching ratio we have used the same parameters as in ||. 
They are given in Table 2. In addition we have included small l/m| corrections as in || and 
also a 3% enhancement [|ItJ from \jvr? c corrections |]][] which were not known at the time of the 
analysis [0]. The relative importance of various uncertainties is shown in Table 3. Comparing 



3 In contrast to a 3% suppression found originally in H (except for the second paper which actually discusses 
the exclusive channels), the corrections to the B ~ * X s j decay have been shown to be positive in 1 10 . 
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Table 2. Input parameter values and their uncertainties. The masses are given in GeV. 



this table with the first row in the corresponding table in || we observe that except for the scale 
uncertainties discussed above, our error analysis agrees well with the one presented in ||. 



Scales 


ol s {M z ) 


m>t,pole 


?Tlc,pole/ TTlbtfole 


TTlb,pole 




CKM angles 


B — > X c ev e 


±4.8% 


±2.9% 


±1.7% 


±5.4% 


±0.7% 


±1.8% 


±2.0% 


±3.8% 



Table 3. Uncertainties in Br(B — > X s y) due to various sources. 



Adding all the uncertainties in quadrature we find 

Br(B^X s j) = (3.60 ± 0.17 (scale) ± 0.28 (par)) x 10" 4 = (3.60 ± 0.33) x 10" 4 (37) 

where we show separately scale and parametric uncertainties. 

Comparing this result with the one of [|J as given in @ we observe that in spite of the smaller 
scale uncertainties in (|37|) our final result is compatible with the one of || and the one given in 
This is due to the parametric uncertainties which dominate the theoretical error at present. Once 
these parametric uncertainties will be reduced in the future the smallness of the scale uncertainties 
achieved through very involved QCD calculations, in particular in []3], |4|, 2C, 21, 11, 12, 13], can 



be better appreciated. This reduction of the theoretical error in the Standard Model prediction 
for Br(B—*X s '~f) could turn out to be very important in the searches for new physics. To this 
end also a better understanding of non-perturbative corrections |5] beyond those considered here 
should be achieved. 

The theoretical estimate in (|37D is somewhat higher than the CLEO result in (|l]) and rather 
close to the ALEPH result in (0). In any case we conclude that within the remaining theoret- 
ical and in particular experimental uncertainties, the Standard Model value is compatible with 
experiment. It will be exciting to watch the improvements in the theoretical estimate and in the 
experimental value in the coming years. 
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